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Abstract 



We study the physics of spontaneously broken family symmetries acting on 
the third generation. Massless familons (or Majorons) / associated with such 
broken symmetries are motivated especially by cosmological scenarios with 
decaying tan neutrinos. We first note that, in marked contrast with the 
case for the first two generations, constraints on third generation familon 
couplings are poor, and are, in fact, non-existent at present in the hadronic 
sector. We derive new bounds from B^-B^ mixing, svv, 
and astrophysics. The resulting constraints on familon decay constants are 
still much weaker than those for the first and second generation. We then 
discuss the promising prospects for significant improvements from searches 
for T^lf, B ^ {T^,K)f, and b {d,s)f with the current CLEO, ARGUS, 
and LEP data. Finally, we note that future constraints from CLEO III and the 
B factories will probe decay constants beyond 10^ GeV, well within regions 
of parameter space favored by proposed scenarios in neutrino cosmology. 
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I. INTRODUCTION 



For over half a century, one of the major puzzles in particle physics has been the question 
of why quark and lepton families replicate. Although we have accumulated a wealth of data 
concerning the masses and mixings of quarks and leptons, we still appear to be far from 
a true understanding of family structure. In the absence of a concrete model to consider, 
it is natural to postulate the existence of some family symmetry that plays a role in 
determining the observed particle spectrum. Once we consider such a family symmetry, we 
face a plethora of options. The symmetry may be (1) discrete,^ (2) continuous and local, or 
(3) continuous and global. Within each of these categories, one may choose any of a number 
of symmetry groups, and the overall family symmetry may even be a combination of the 
three possibilities. 

Of course, any exact family symmetry of the underlying theory must be spontaneously 
broken at some energy scale since we know that the quark and lepton masses are very 
different from one family to the next. For option (1), spontaneously broken discrete sym- 
metries, domain walls are the only model-independent predictions, and these cannot be 
studied in particle physics laboratories. In case (2), the masses of the family gauge bosons 
of spontaneously broken local continuous symmetries can be constrained, e.g., from K^-K^ 
mixing 

From a phenomenological point of view, however, possibility (3) is particularly enticing, 
as it implies the existence of massless Nambu-Goldstone bosons, called "familons," from the 
spontaneously broken family symmetry. This family symmetry may be either Abelian or 
non-Abelian; Nambu-Goldstone bosons associated with the spontaneous symmetry breaking 
of an Abelian lepton number symmetry are often called "Majorons."0 The existence of new 
massless particles has many implications in particle physics, astrophysics, and cosmology, 
and, as we will see, may be probed in a wide variety of experiments. Moreover, the cou- 
plings of familons at low energies are determined by the non-linear realization of the family 
symmetry. These couplings are, e.g., of the form 

jd,r , (1.1) 

where F is the family symmetry breaking scale, i.e., the familon decay constant, /° are 
the familons, are the generators of the broken symmetry, and the ipi are fermion fields 
in terms of which the flavor symmetry is defined. The strength of the familon coupling is 
therefore inversely proportional to F and can be constrained for a given family symmetry 
group in a model-independent manner. 



^ There is a subtle distinction between global and gauged discrete symmetries For this phe- 
nomenological analysis, however, they are equivalent. 

^Majorons have been extensively studied, and arise in a variety of models including, for 
example, supersymmetric theories with spontaneous -R-parity breaking In this paper, we study 
a number of probes, many of which are applicable to both Abelian and non-Abelian symmetries. 
We use the generic name "familon" to denote the associated Nambu-Goldstone bosons in either 
case. 
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Familon couplings between the first and second generations have been studied extensively 
and will be reviewed below. In contrast, however, couplings involving the third generation are 
largely unexplored, although they may have rather rich phenomenological and cosmological 
implications P]. Current constraints in the lepton sector are relatively weak, with the best 
bounds coming from r —>■ {e,fi)f bounds 0, and there are at present no corresponding 
bounds reported in the hadronic sector (see, however, Ref. At the same time, it is a 

logical possibility that the familon couples preferentially to the third generation, and models 
have been proposed in which this is the case |ll|]. It is therefore interesting to explore the 
possibilities for improving (or setting) bounds on familon scales for the third generation, 
especially in light of the upcoming B physics experiments. 

In this paper, we will study what we believe to be the most sensitive probes of couplings 
of familons to the third generation, primarily to r leptons and b quarks. We show that 
dedicated analyses of existent data from CLEO, ARGUS, and LEP could probe family 
symmetry breaking scales up to ~ 10^ GeV and may be significantly improved at future 
B factories. Simply because this is largely unexplored physics, there is a high discovery 
potential for familons at these facilities. 

Familon couplings to the third generation are also of interest from a cosmological point 
of view. The mass of the r neutrino is still allowed to be as large as 18.2 MeV experimen- 
tally |jl2[ . A heavy r neutrino has interesting consequences for both big-bang nucleosynthesis 
(BBN) [|13|-|1^ and large scale structure formation |ll7H20|], as will be discussed in Sec. 0. 



Since a heavy neutrino (> lOO/i^ eV, where h is the expansion rate of the universe in units 
of 100 km/sec/Mpc) must decay in order not to overdose the universe, an invisible decay 
into a lighter neutrino and a massless boson, such as a familon (or Majoron), is typically 
required. (The three neutrino mode is strongly disfavored and therefore the familon mode 



is most preferred |T3[.) There is therefore an interesting interplay between experimental 
searches for familons and scenarios requiring heavy neutrinos, and, as we will see, future 
collider experiments and analyses may severely constrain a number of such cosmologically 
motivated scenarios. 

This paper is organized as follows. We begin in Sec. |T| with a discussion of familon 
interactions. In particular, we emphasize that the familon interactions of particles in the 
same gauge multiplet are expected to be comparable. In Sec. |T| we consider constraints 
on familon interactions that may be inferred from current experimental data, concentrating 
on familon couplings to the third generation. Current bounds on third generation couplings 
from astrophysical considerations are presented in Sec. |rV|. We then describe some promising 
prospects for detecting familons in B physics at future experiments in Sec. 0. Finally, we 
note some of the interesting cosmological implications in Sec. |V| and give our conclusions 



in Sec. VII 



II. FAMILON INTERACTIONS 



The standard model contains 15 particle states in each of the 3 generations. These 
states are distinguished by the SU(3)cxSU(2)i,xU(l)y gauge interactions, which divide 
each generation into 5 multiplets: Q, U, D, L, and E. The gauge interactions therefore break 
the flavor symmetry group from U(45) to U(3)^. In the standard model, the flavor group 
U(3)^ is broken explicitly to U{1)b x U(l)i by Yukawa couplings. However, in extensions 
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of the standard model in which one hopes to gain some understanding of the pattern of 
fermion masses and mixings, some subgroup of the flavor group may be an exact symmetry 
of the Lagrangian that is broken spontaneously by the vacuum, and it is this possibility we 
consider here. 

The massless Nambu-Goldstone bosons of the spontaneously broken flavor symmetry, 
familons have interactions given by the couplings 

Cf = ^d,rj^\ (2.1) 

where are the familon fields, and are flavor currents.^ The interactions are suppressed 
by F, the scale at which the flavor symmetry is spontaneously broken. Note that familons 
are derivatively coupled,0and so do not mediate long-range (~ r~^) forces. The most general 
current J'^"' composed of two fermion fields takes either the form 

J'''' = i'a''i9v + 9A7,)Ttj^l'j (2.2) 

or 

= ^Pa^igLPL + gRPR)m^j , (2.3) 

where Pl,r are the projection operators |(1 ± 75), i and j are generational indices, and T^" 
are the spontaneously broken generators of the family symmetry. The fields ipi and i/jj are 
fermion mass eigenstates, which we assume here to be also flavor eigenstates. (The more 
general case is described below.) Using the form of the current given in Eq. (|2.2| ), the familon 
interaction may be written as 

= -^Pi'i [gv{fni - TUj) + gAirrii + 772^)75] ^ipj , (2.4) 

where in the last step we have integrated by parts and then substituted the equations of 
motion. The second line of Eq. (|2.4| ) is of course only valid for on-shell fermions such 
as external leptons, whereas in hadronic matrix elements and processes including off-shell 
fermions, the derivative coupling of the first line must be used. 

We see that familons may mediate or be produced in family-changing processes. They 
may also couple to identical fermions ipi = ipj, but only through axial couplings. What 
processes are mediated by familons depends on the particular family symmetry group that 
is broken. For example, for O(A^) groups, the generators Tjj are anti-symmetric, and so 
do not generate flavor-diagonal interactions. However, they do generate interactions like 



^Throughout this study, we will assume that no additional light degrees of freedom are introduced 
by other new physics. 

^If the flavor symmetry is anomalous, familons may also have non-derivative, flavor-diagonal 
couplings. We will not consider such couplings here. 
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f'4'il54'j~f'4'jl5'^iy wheie we have considered axial vector current interactions as an example. 
Familons from O(A^) groups may therefore mediate neutral meson mixing, which we will 
consider in Sec. |111 B\ The situation is reversed for SU(iV) groups. Here, flavor-diagonal 
couplings exist. However, if we consider any SU(2) subgroup and form the complex familon 
f = +ip, the off-diagonal interactions are given by fipi'j^ipj + f*i^j'~i^ipi, and we see that 
/ exchange cannot induce neutral meson mixing. 

Up to this point, we have ignored possible mass mixing effects. In general, if the flavor 
eigenstates ip' are related to the mass eigenstates ip by 

i^' = U^^ , (2.5) 

where [/^ is a 3 x 3 unitary mixing matrix, the familon interactions are given by 

= jd,r ^^l^{9v + 9Al^)Th'^, , (2.6) 

where = U\T'^U^. Mass mixings may therefore generate flavor- diagonal interactions 
from flavor off-diagonal interactions, and vice versa. For example, in the case of an Abelian 
U(l) symmetry, mass mixing effects may generate flavor- changing interactions. They may 
also extend non-maximal family symmetries to couplings involving all three generations; for 
example, a U(2) symmetry between the first and second families, may, after rotation to mass 
eigenstates, result in familon interactions involving the third generation. 

While the phenomenology of familons varies from group to group, it is important to 
note that gauge symmetry relates the familon interactions of particles in the same gauge 
multiplet. As an example, let us consider a spontaneously broken lepton flavor symmetry. 
The familon interaction is then given by 

^^d.rl'a^Tt.L] , (2.7) 

where the SU(2) lepton doublets L\ = (z/^', are in the flavor eigenstate basis. This interac- 
tion therefore generates familon interactions for both the charged leptons and neutrinos. In 
the presence of neutrino masses, the flavor eigenstates may not correspond to mass eigen- 
states. The familon interactions in the mass basis are then 

^fd,r ^.7^t;,^-^, + ^fd,r ui'n.h , (2.8) 

where = UlT^U^, and we have defined V = U^Ui and T;"^ = V^T;^V. and are there- 
fore related by a similarity transformation, and in the presence of mass mixing, the couplings 
of the interactions of d^f^ i^n'^i^j and 9^/" h'j^lj are not necessarily identical. However, in 
the absence of fine-tuning, we expect these couplings to be of the same magnitude. Bounds 
on one familon interaction may thus be considered to imply comparable bounds on the other 
interactions linked by gauge symmetry. 

Because the familon interactions of particles in the same gauge multiplet are comparable 
in the absence of fine-tuning, there are many more relations in theories with enlarged gauge 
groups. For example, for SU(5) grand unified theories (GUTs), the particles dji, v^cl C 5 are 
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expected to have comparable familon interactions, as are the particles ULidi^UR^eR C 10. 
A particularly relevant example for our study below is that, in the GUT framework, bounds 
on familon decays of B mesons imply bounds on familon decays of tau neutrinos in the 
absence of fine-tuning. 

Flavor mixing effects also induce familon couplings of fields with different generational 
indices. In the quark sector, for example, substituting the quark doublet Q[ = (m^, d'^ for L[ 
in the discussion above, Eq. (|2.8D becomes 



^d,r ua'n^^u, + ^;d,r da'T^.,d, , (2.9) 

where and are related by the CKM matrix through 

= ViKMT:VcKM . (2.10) 

We see that in general, couplings to all generations are induced by flavor mixings. For 
example, a familon with flavor-diagonal coupling to tt in the up sector couples not only to 
bb, but also to, for example, bs and dd. The induced couplings to first and second generation 
quarks in this case are CKM-suppressed, but may still lead to significant bounds when, as 
is often the case, these induced couplings are much more strongly constrained. We will 
consider the constraints on mixing- induced couplings from K decays in Sec. [Ill A| and from 



supernova cooling in Sec. |1V C 



Finally, note in Eq. (p.2|) that the strength of the interaction depends not only on F, 
but also on T^j and the couplings gv,A- In the following sections, we will present a vari- 
ety of bounds on combinations of these couplings, and it is important that we define our 
conventions and normalizations. We will always define our interaction as 

jdJ4j,Y{9'^+9h5)^Pj , (2.11) 

and similarly for g^l and (7^; the superscripts of the couplings will often be omitted when 
they are obvious from the context. In presenting our bounds, it will be convenient to define 

F!,^F/gy, (2.12) 

where I = V,A, L, R. In addition, as many of our bounds are to a good approximation 
independent of the chirality of the interaction and so only dependent on the combination 
(7y^ + we define 

Fij = , ^ . (2.13) 

^ ij2 , ij2 ^ ' 

'Qv +9a 



III. BOUNDS FROM ACCELERATOR DATA 

As described in the previous section, familons may take part in flavor- changing pro- 
cesses, and bounds on such processes lead to lower bounds on the familon energy scale. 
For familons mediating transitions between the first and second generation, such bounds 
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are rather stringent. In contrast, similar bounds involving the third generation are much 
weaker, with the previously reported constraints limited only to bounds from rare r decays. 
We are thus motivated to focus on the third generation. In Sec. [Ill A| , we begin by reviewing 
and contrasting such bounds, and then discuss the implications of flavor eigenstate mix- 
ings. We then go on to derive new bounds from a variety of processes. In Sec. |III B| we 
consider familon-mediated processes such as neutral meson mixing and rare leptonic decays 
of mesons. Finally, in Sec. |111 C| we consider possible analyses at LEP and extrapolate a 
preliminary ALEPH bound on 6 ^ sui? to a bound on 6 — > sf. 



A. Decays to familons 

We begin by considering bounds from decays of mesons and leptons to familons. Nor- 
malizing the relevant familon scale according to Eq. (|2.11|) , we find 



nK--^f) = T^'^9'v(3'mO)\\ (3.1) 

where /3 = 1 — m^/m\. In the limit of exact flavor SU(3) symmetry, the form factor 
{tx^ {j)')\s'^^ d\K^ {p)) = Fi{q^){p+p')^ at zero momentum transfer has a fixed normalization, 
-Fi(O) = 1. For leptonic decays li —>■ Ijf, the exact tree-level partial decay width in the limit 
of massless Ij is given by 

m - IJf) = ^^(9'v + 9\)0' , (3.2) 

where here (3 = 1 — mf, / mf_ . 

The strongest bound on any flavor scale is derived from the constraint on exotic K 
decay. Using the above expressions, the experimental result B{K~^ tt^/) < 3.0 x 
10"^° (90% CL) [H leads to the bound 

Fj^ > 3.4 X 10" GeV . (3.3) 

Note that the limit on B{K^ tt^/) bounds only the vectorial familon coupling; the axial 
coupling is unconstrained. For the leptonic sector, Jodidio et al. report the constraint 
B{ii^ —> e+/) < 2.6 X 10"^ (90% CL) which they obtain under the assumption of a 
vector-like familon coupling. This can be converted into the bound 

> 5.5 X 10^ GeV . (3.4) 

For familon interactions of arbitrary chirality, the slightly weaker constraint 

F^e > 3.1 X 10*^ GeV (3.5) 



may be obtained from the bound B{fi^ e'^lf) < 1.1 x 10"^ (90% CL) 

We now compare these bounds to those available in the third generation. The ARGUS 
collaboration has bounded the branching fractions of r decays into light bosons and found 
the limits B{t' fi~ f) < 4.6 x lO'^ (95% CL) and B{t- f) < 2.6 x 10^=^ (95% CL). 
These imply the following constraints on the flavor scale: 
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Frf, > 3.2 X 10^ GeV (3.6) 
Fre > 4.4 X 10^ GeV . (3.7) 

We see that the bounds on flavor scales in the leptonic sector are significantly less strin- 
gent for third generation couplings than for those involving only the first two. The discrep- 
ancy is even more pronounced in the hadronic sector, where there are as yet no reported 
bounds on flavor scales from B decays. 

It is also worth noting, however, that strong bounds on a particular flavor scale, such as 
the one on Fj^, may imply significant bounds on other flavor scales as well. These bounds are 
induced by the flavor-mixing effects discussed in Sec. || and are thus mo del- dependent. As 
an example let us now assume that flavor and mass eigenstates coincide for up-type quarks. 
A given familon coupling in the up sector requires, by gauge invariance, a corresponding 
coupling in the down sector. For example, from Eqs. ( |2.9| ) and ( ^.10[ ) we see that the coupling 



d^f I'-^^Plc/ Fl"^ induces the coupling V*gVcdd^f S'j^PLd/ F^^, which mediates the rare decay 
rr^f. Assuming complex familons, the Hermitian conjugate coupling gives a similar 
contribution oc V*gVtd to the decay into the complex conjugate familon. Summing both 
decay widths and comparing to the bound on F^^ in Eq. ( p.3|) , one can derive the mixing 
induced bound 

> 2.2 X 10^ GeV . (3.8) 
Under similar assumptions, we find F^^ > 6.6 x 10^ GeV. Note, however, that such bounds 



do not apply if the mass and flavor bases are aligned in the down sector pj] or if the 
couplings are purely axial. 

B. Familon-mediated processes 

In this section we derive new constraints on the scale of spontaneous flavor symmetry 
breaking by considering non-standard familon contributions to neutral meson mixing and 
existing bounds on rare leptonic decays such as re. 

A familon contribution to neutral meson mixing requires a real flavor group to be spon- 
taneously broken in the corresponding sector, such that the same real familon scalar field 
couples to the quark current and its Hermitian conjugate. For concreteness, let us consider 
the B^ — B^ system; similar formulae hold (at least approximately) for other neutral meson 
systems. Assuming the general coupling structure 

jdj[dr{9v + 9Al5)b-bY{gv + 9Al5)d\ , (3.9) 
we find a familon contribution to the mass splitting of 

Amif) ^ \m,o - mso\ ^ ^Sl^ . (3.10) 
Eq. ( |3.1(J| ) may be derived by taking the matrix element of the non-local operator 

1 1 — — 

^J^d^l^'iav + gAl^)hi3Afd^Y{9v + QAlblhs (3.11) 
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TABLE I. Bounds on the flavor scale from contributions to neutral meson mixing from familon 
exchange as given in Eq. ( |3.1[1| ). Note that these limits do not apply to vector-like couplings, and 
that this process requires a real flavor group so that a real familon scalar fleld couples to a current 



operator and its Hermitian conjugate, as in Eq. (3.9). 



Am, 



'cxp 



Bound 



0.5 X 10^2 ^j^-i 

< 21 X 10^° hs-^ |7| 
0.53 X 10^^ hs-^ |5| 



F^^ > 6.4 X 10^ GeV 
> 6.9 X 10^ GeV 
F^^ > 1.7 X 10^ GeV 



between and states and using the definition of tlie pseudoscalar decay constant, 
(0|67^75(i(0)|-B°(p)) = if^op^. Tlie subscripts a, /?, 7, and 6 in Eq. are color indices. 

Between two color singlet states, tliere are two contributions. Tlie first one arises from 
a = f3 and 7 = 5 witli a familon in tlie s-cliannel. In this case, the momentum transfer 
through the familon propagator is = m^o , and after a vacuum insertion, it is easy to verify 
that this contribution is as in Eq. (|3.10|) , but without the factor of 5/6. However, there is 
also a t-channel contribution from a = 6 and /5 = 7, which may be evaluated by a Fierz 
transformation and then a vacuum insertion as before. For a heavy-light system like the i?° 
meson, one may assume the free-quark picture, in which the momentum transfer is governed 
by the energy of the "static" b quark g° irib ~ m^o, and, in the numerator, the derivative 
acting on the quark current gives again a factor of mi,. Using {0\b'y^d{0)\B'^{p)) ^ if^orriBO 
and including the relative color factor of 1/3, one can estimate the t-channel contribution 
to be —1/6 times the s-channel contribution, which leads to Eq. ( |3.1(J| ). 

Our result should be fairly reliable for the B^ meson. For and mesons, the 
evaluation of the t-channel momentum transfer is more ambiguous. However, because this 
contribution is suppressed relative to the s-channel part, we expect the result of Eq. ( p.lO| ) 
to be reasonably accurate in these cases as well. We also note that a vector-like familon 
interaction does not contribute to the mass splitting, at least in the heavy quark approx- 
imation rrib ~ rriBO. Although one might expect a vector contribution to appear in the 
t-channel contribution after the Fierz rearrangement, one finds that the term proportional 
to gv contains axial vector and pseudoscalar contributions of equal magnitude but opposite 
sign. 

The constraint on the flavor scale F results in principle from the requirement that the 
combined standard model and familon contributions do not exceed the measured value. 
However, when considering nonstandard contributions, it is also uncertain what one should 
take as the standard model contribution. For example, the reported value p5[ for iV^^^dl 
is derived from B^ — B^ mixing under the assumption that the standard model gives the 
only contribution. As a conservative bound, we simply compare the familon contributions 
directly to the corresponding measured values. The results are summarized in Table |. We 
take the decay constants to be fso ~ 175 MeV and /^o f» 205 MeV from recent lattice 
results and f^o ~ fK+ ~ 160 MeV [^. Since we use the measured mass splitting (not 
its error), the bounds from B^ and will only improve when the size of the standard model 
contribution can be quantified independently. For the Z)°, where only the upper bound on 
the mass splitting is known, future experiments will improve the bound. 

We next consider rare leptonic decays of neutral mesons, mediated by familon exchange. 
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TABLE II. Limits on flavor scales and couplings for some rare meson decays. The branching 
ratio bounds are on the sum of the two charge states, assuming real familon scalars that mediate 
both decay modes. If familons mediate only one decay mode, the quoted bounds on F are weakened 
by a factor of 2^/^. In calculating these bounds, we neglect small corrections from the lighter lepton 
mass. 



Branching Ratio Upper Bound 



Bound 



5.3 X 10-"^ H 



{F^^FreYI"^ > 3.5 X 10^ GeV 
(F,>,^)i/2 > 3.1 X 103 GeV 
(Fft^F^e)'/' > 2.8 X 103 GeV 
{FiF^,)^^ > 1.2 X 103 GeV 
^ii^/.e)^/^ > 3.8 X 10^ GeV 



8.3 X 10"^ Us] 
5.9 X 10-6 l28| 
1.9 X 10-5 |29| 
3.3 X 10-11 |C| 



Such decays are possible if the same familon couples to both quarks and leptons. This 
is guaranteed in grand unified scenarios, where quarks and leptons are in the same gauge 
multiplet. In general the relevant interaction can be written in terms of effective vector 
and axial vector couplings that parametrize the familon couplings and mixing angles of a 



particular model. 
Lagrangian 



For example, the process B — > r+e can be mediated by the interaction 



1 



+ h.c. 



(3.12) 



Note that the constants gv and gA may be different in the hadronic and leptonic sectors. 
Also, even if familon couplings always include third generation flavor eigenstates, mixing 
effects may induce transitions like fi'^e~ . 

With the interaction defined in Eq. ( p.l2| ) one obtains a width of 



T e 



1 f^og^^^mBom: 



{9r+9r)f3'-2'^{gr-9r 



(3.13) 



where (3 = 1 — m^/m^o, and we have displayed the leading gy — g\ piece. In the limit where 
the lighter lepton is massless, the result is independent of the chirality of the interaction and 
depends only on the combination of lepton couplings gy + g\- Expressions for other similar 
processes are obtained by replacing the coupling constants gy, gA and the meson and lepton 
masses accordingly. Limits on the flavor scales from current experimental bounds on rare 
leptonic decays are given in Table 0. 

The bounds of Tables | and |I| are significantly weaker than those presented in Sec. [Ill A . 
This is especially true in Table |I[ as rare leptonic meson decays are dependent on the flavor 
scale to the fourth power. However, such processes set bounds on third generation hadronic 
familon couplings, which were previously unconstrained. It is also important to note that 
the bounds on familon couplings to the first two generations are also interesting, as they 
constrain axial couplings, whereas the bound from K decay reviewed in the previous section 
bounds only vector-like couplings. 
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C. Constraints from LEP 



Currently, there are no reported experimental bounds on decays b —>■ {s,d)f. One can, 
however, infer a constraint from ALEPH's preliminary bound on 6 — suu |^T[. By searching 
for events with large missing energy, they placed the constraint B{b —>■ svv) < 7.7 x 10~^ 
(90% CL). One can rescale this constraint to obtain an upper bound on B{h sf). 



The analysis for b svv relies on the -Emiss distribution |32], where -Emiss is defined by 



in each hemisphere of 6-tagged events. Here, -Ebeam is half of the center-of-momentum energy, 
-Ecorr = (M^ame ~ M^pp ) /^-E'beam , where Mgame and Mopp are the visible invariant masses in 
the same and opposite hemispheres, respectively, and i?vis is the total visible energy in the 
hemisphere. -Ecorr improves the estimate of the actual missing energy in the hemisphere by 
correcting for the fact that the hemisphere with larger invariant mass typically has higher 
energy. 

The backgrounds from b — >■ IvX and c — > IvX are suppressed by rejecting events with 
identified or /x^ in the relevant hemisphere. Up to this point, we do not expect significant 
differences in efficiencies between the b svv mode and the b —>■ sf mode. They then 
required 35 GeV < -Emiss < 45 GeV. The efficiencies for this requirement obviously differ 
between the two decay modes, since the mode b — > suu has two missing neutrinos, resulting 
in a harder -Emiss spectrum than that of the b ^ sf mode. The -Emiss spectrum of both modes 
may be calculated by convoluting the theoretical missing energy distribution in three-body 
(suu) and two-body (sf) decays with the measured b fragmentation function [Q. We find 



that the ratio of efficiencies is 0.43 with little dependence on the details of the fragmentation 
function. By scaling the reported B{b —>■ svv) upper bound by this factor, we find 

5(6 ^ s/) < 1.8 X 10"^ . (3.15) 



Using the expression of Eq. ( |3.2| ) with the substitution of m^o mi, for m/^, this corresponds 
to a limit on the flavor scale of 

Fb, > 6.1 X 10' GeV . (3.16) 

Note that this analysis does not require an energetic strange particle, and so the constraint 
of Eq. (|3.15|) is actually on the sum B{b s f)-\-B{b ^ d f). Thus, for ~ F^, the bound 



on the flavor scale given in Eq. ( |3.16| ) improves by a factor \/2. The bound of Eq. ( |3.16|) is 
enhanced by the fact that the SM decay width is greatly suppressed by V^b, which increases 
the sensitivity of b decays to small exotic decay widths. 



IV. BOUNDS FROM ASTROPHYSICS 

In this section, we discuss constraints on third generation familon couplings from astro- 
physics. We begin in Sec. |1V A| with constraints on direct (tree-level) couplings. Second 



and third generation particles are absent in almost all astrophysical objects. The exception 
is supernovae, where all three neutrino species are thermalized in the core. We therefore 
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consider what bounds on familon couplings to r neutrinos may be obtained by supernova 
observations. Couplings of familons to the third generation may also radiatively induce cou- 
plings to first generation particles. Although such induced couplings are suppressed by loop 
factors, they are so stringently bounded by constraints from supernovae, white dwarfs, and 
red giants that interesting bounds also result. These are studied in Sec. |1V B[ Finally, mix- 
ings of flavor eigenstates may also induce couplings of familons to the first generation; such 
effects are discussed in Sec. [IV C| . It is important to note that, while the bounds derived in 
this section are rather strong in certain cases, they are also typically more model-dependent 
than, for example, the accelerator bounds of the previous section. We therefore specify the 
necessary conditions for each bound in detail in each case. 



A. Bounds from direct couplings 

In 1987, the Kamiokande group and the 1MB group independently detected neutrinos 
emitted from supernova SN 1987A. They observed that the neutrino pulse lasted for a few 
seconds. Furthermore, their results indicate that neutrinos carried off about 10^^ erg from 
the supernova. The observed duration time and neutrino flux can be well explained by the 
generally accepted theory of core collapse, and the observations confirmed the idea that most 
of the released energy in the cooling process is carried off by neutrinos. Exotic light particles, 
such as familons, may affect the agreement of theory and observation, since they can also 
carry off a significant energy fraction. The core of the supernova is hot (T ~ 30 MeV) 
and dense, and so neutrinos are thermalized in the core and can be a source of familon 
emission. If the energy fraction carried away by familons is substantial, the duration time 
of the neutrino pulse becomes much shorter than the observed value. In order not to affect 
the standard cooling process, the familon luminosity Qj must be smaller than the neutrino 
luminosity, i.e., less than ~ 10^^ erg/sec. 

This constraint can be satisfied in two different regimes of the familon coupling strength. 
For sufficiently high flavor scales F, the familon interaction is weak enough that familons 
are rarely produced and the familon luminosity Q/ is suppressed. On the other hand, 
for sufficiently low flavor scales, although familons are readily produced, they interact so 
strongly that they become thermalized and trapped in the core as well, thus decreasing the 
familon luminosity. Therefore, there are two parameter regions consistent with observations, 
high and low F, with an excluded region in the middle. 



These types of constraints have been discussed by Choi and Santamaria |33] in the con- 
text of a Majoron model. We modify their discussions slightly for the familon case. To 
simplify the analysis, we will look at two extreme scenarios. First we consider a diagonal 
familon coupling to v^, as in the case of an Abelian family symmetry, and second we analyze 
a purely off-diagonal coupling, as in the case of an 0(2) family symmetry.^ For a general 



^Throughout our discussions, we assume that neutrinos are Majorana particles. Observations of 
supernova SN1987A imply that Dirac neutrinos must be Ughter than 3 keV |36| or heavier than 
31 MeV As we will discuss in Sec. VI, most of the interesting mass range from a cosmological 
point of view is therefore excluded. 
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family symmetry, one expects familons with both diagonal and off-diagonal couplings; a 
generalization to such cases is straightforward. In this subsection, we also neglect possible 
mismatches between the flavor and mass eigenstates, and assume that the relative angles re- 
lating the two are small. Such mismatches will be discussed in Sec. [IV Q . Finally, we assume 
that rriu^, m^^ are negligible compared to m^^, as suggested from laboratory constraints as 
well as the corresponding masses of the charged leptons. 



Familon with diagonal coupling 



Here we consider a purely diagonal familon coupling to v-r-, such as in models with a 
U(l) family symmetry acting on the third-generation lepton doublet {vtiTl}- The relevant 
interaction is given by 



(4.1) 



Let us first consider the case where the familon can freely escape the core of the supernova. 
Based on the interaction given in Eq. ( f4.1| ), potentially significant processes of familon 
production are the neutrino scatterings Vr^r f f and Ur Vrf , the latter process being 
allowed due to background matter effects. The familon luminosities due to these processes 



are given in Ref. [p^ : 

Qf{vri^^ ff) ^ 8.8 X 10^3 erg/sec x 



VMeV 



' \^ VrVr J 



1.6 X 10^^ erg/sec x 
1.9 X 10^° erg/sec x 



'MeV\ VGeV\^ 



/GeVV 



m, 



, m^^ > 95 keV 



(4.2) 



(4.3) 



, m^, < 95 keV 



If either one of the above luminosities is larger than ^ 10^^ erg/sec, the cooling process of 
the supernova may be dominated by familon emission, and the duration time of the neutrino 
pulse becomes shorter than 0(1 sec). Imposing the constraint that the familon luminosities 
given in Eqs. ([4.2|) and (|4.3|) are smaller than 3 x 10^^ erg/sec, we obtain the following 



constraints: 



ff-- 



^rf ■■ < 



/GeV\' 



'MeVy /GeV\ 



< 5.8 X 10"^ 



< 0.43 



m^^ > 95 keV 



m 



(4.4) 



(4.5) 



Familon volume emission is sufficiently small when both Eqs. ( [4.4| ) and ( [4.5[ ) are satisfied. 

On the other hand, if familon interactions are strong enough, familons effectively scatter 
off the neutrinos in the background and get thermalized and trapped in the supernova. Once 
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this happens, a thermal sphere of familons is formed, just hke the thermal neutrino sphere, 
and familons can only be emitted from the surface. The familon luminosity essentially obeys 
the formula of blackbody emission with the surface temperature of the familon sphere. The 
important point is that, once the familon is trapped, the familon luminosity decreases as the 
familon interaction becomes stronger. This can be understood in the following way: as the 
familon interaction becomes stronger, familons can be thermalized with a lower temperature. 
(Notice that the scattering rate increases for higher temperature.) The surface temperature 
of the familon sphere then decreases, and hence the luminosity is suppressed. Therefore, 
the familon luminosity can be small enough when the scale F is sufficiently low. Following 
Ref. we find that the cooling through familon emission is sufficiently suppressed {i.e., is 
less than 3 x 10^^ erg/sec) when either one of the following constraints is satisfied: 



f \ fGeVV 

-4; > 8.3 X 10 

r_\(GeV 



V. \ I'tI't / 



(4.6) 



(4.7) 



Of course, if the familon has strong interactions with other light particles (the photon, 
electron, or neutron), familons may be trapped by other processes as well. This gives 
additional regimes where the earlier constraints of Eqs. ( [4. 4] ) and ( [4.5| ) can be evaded. 

In Fig. |l|, we show the upper and lower bounds on the diagonal coupling -Fjj;^^ as a func- 
tion of the neutrino mass m,^^. The region above the upper line is allowed because familon 
emission is sufficiently suppressed by the flavor scale F. This line is basically determined 
by Eq. ( [4.4|) ; the slight bump is due to Eq. ( [4.5|) . As we can see, the lower bound on F is 
at most 1 TeV for the maximum allowed value of the tau neutrino mass (18.2 MeV), and it 
becomes less stringent as the mass becomes smaller. The lower boundary is determined by 
Eq. (|4.(j|), which supercedes Eq. (|4.7|). The region below this line is also allowed because the 
familon is trapped in the core and the contribution to the cooling is again sufficiently small. 



Familon with off- diagonal coupling 

Here we discuss SN 1987A constraints on a familon which has only an off-diagonal cou- 
pling, such as in the case of an 0(2) family symmetry. The low-energy Lagrangian of the 
model can be written as 

Cf = jgr^d.f {vrrPLy, - y.YPLyr) ■ (4.8) 

The inclusion of Ue in the discussion is straightforward. The familons are produced by 
VrVr ff via t-channel v^ exchange, v^Vf^ ff via t-channel Vr exchange, or the decays 
i^T ^tif ■ Following Ref. again. 



QAi^ri^r - //) - 2.2 X 10- erg/sec x (^)' , 



(4.9) 
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FIG. 1. Excluded region in Fj;'^^^ as a function of the neutrino mass rrii,^, derived from SN1987A. 
The bounds shown correspond to the case of diagonal familon coupling. For the off-diagonal case 
the bounds on F^;'^^^ are very similar: the excluded region is only marginally shifted down by a 
factor of 1.2, and the little bump from v-j- v-j-f is absent. 

Qf{y,v, ff) ^ 2.2 X 10*^=^ erg/sec x (^)' (^^) , (4.10) 
QA^r - - 3.0 X 10^^ erg/sec x (^)' . (4.11) 

We require that all of these familon luminosities are smaller than 3 x 10^^ erg/sec, and obtain 
the following constraints: 



(£^)(|g)% 8,3X10-. 
fe)fM!X > 3.3X10-=. 



(4.12) 
(4.13) 



where the lifetime of u^- is given by 



1 rnt^ 
167rF,f ,2 



(4.14) 



Increasing the interaction strength further into the excluded region, familons eventually 
become trapped and rendered harmless again. This occurs when any one of the following 
constraints are satisfied: 
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i";^^"^' ■■ f^) fl^V > X 10- , (4.15) 

fi^, ^ i^r : {r.Jsec) (^] < 10"^ . (4.16) 

As mentioned before, if the familon has strong interactions with other hght particles, these 
interactions may lead to thermalization of familons as well, resulting in additional allowed 
regions for low flavor scales. 

The resulting excluded region is fairly similar to that of the diagonal coupling case. 
The constraint from the decay process Eq. ( [4.13| ) is important only for smaller F^^^, or 



larger masses my^^ values that are outside our range of interest. In addition, the small 
bump in Fig. |l] now disappears due to the absence of the v.^ Urf process. The dominant 
constraints are therefore from Eqs. ( [4. 121) and ( [4.15|) in the off-diagonal case, which differ 



from the dominant constraints of Eqs. ( [4.4|) and ( |4.6| ) in the diagonal coupling case only by 
a small constant factor. The boundary of the excluded region is therefore given by the lines 
of Fig. |1| shifted downwards by a factor of 1.2 in F. 



B. Bounds from loop- induced couplings 

In the previous subsection, we considered astrophysical constraints on tree-level familon 
couplings to Ut-. In addition, however, astrophysical bounds may also be used to constrain 
familon couplings to all other particles, as these couplings may induce couplings of familons 
to electrons and nucleons at the loop level. While these induced couplings are suppressed 
by the usual loop factors, the bounds on familon couplings to flrst generation particles are 
so stringent that these constraints may be strong in certain cases. In fact, we will see 
below that the contributions to induced couplings are proportional to fermion masses, and 
so these constraints are particularly relevant for couplings of familons to third generation 
fermions. In this subsection, we will estimate the induced couplings for various choices 
of the family symmetry group and determine what lower bounds on flavor scales F result 
from current astrophysical constraints. For simplicity, we will limit our discussion here to 
familons with flavor- diagonal couplings to the third generation, and ignore possible rotations 
relating the flavor and mass eigenstates. Extensions of this analysis to more general cases 
are straightforward. 

To evaluate the strength of the induced coupling, we will begin by considering the low 
energy effective theory below the flavor scale F. In this approach, the theory is specifled 
by the flavor symmetry, that is, the low energy derivative couplings of the familon, and no 
further knowledge of the mechanisms of flavor symmetry breaking is required. With the 
assumptions given above, the dominant contribution to the induced couplings is from the 
Z-f mixing graph shown in Fig. |^. Here x is any one of the third generation particles 
directly coupled to the familon, and ip = u, or d. (There are also additional contributions 
from penguin-like W diagrams, but these are suppressed by mixing angles, e.g., V^^ in the 
case oi ip = d.) Let us deflne the f-x coupling as 

jd^fxr{9LPL + gRPR)x, (4.17) 
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FIG. 2. The Z-j mixing diagram. 



and the Z-x coupling as 



{91Pl + 91Pr)x 



(4.18) 



where gl = gz^I^ - Q^, 



sm' 



9l 



—gzQx^^^'^^w, and gz = e/(sin6'vi/ cos6'iy). The 



induced Z-f mixing from the fermion loop is divergent, and the logarithmically-enhanced 
contribution is 



2Nr 



z-f 



(47r)2F 



In At 



ml{g 



L - gnM - gl)g'"' - \{gLgl + gRgDiv^g'" - p'f 



(47r)2F ^ 



A2 



mi 



{gL-gR){gl-gl)z^d,f 



(4.19) 



where A''^ is the number of colors of the fermion x, and A is the effective ultraviolet cutoff of 
the order of the flavor scale F. Note that gi ~ gn = gz^x fermion charge drops 

out: the Ward-Takahashi identity guarantees that the sin^ 9w piece in the Z vertex gives 
a completely transverse vacuum polarization amplitude proportional to p^g^'^ — p^p", which 
vanishes when contracted with 9^/. The leading contributions to the induced couplings are 
determined by the amount of current non- conservation, i.e., the masses of the particles in the 
loop, and the third generation couplings therefore give the most important contributions.^] 



The mixing of Eq. ( [4.19| ) is logarithmically-enhanced and so typically gives the leading 
contribution if present. However, there are cases in which this term is not present or is highly 
suppressed. First, it may be that the amount of current non- conservation is itself suppressed 
by inverse powers of the flavor scale. For instance, in the singlet Majoron model [^, lepton 
number conservation in the low-energy theory is violated only by neutrino masses which are 
of order 1/ F due to the seesaw mechanism. The neutrino loop contribution to the Z-Majoron 
mixing is then 1/ F^ and highly suppressed. Second, if the familon coupling is vector-like so 
that gL = gRi the logarithmically-enhanced term is absent. For example, a familon coupled 
to the (possibly generation-dependent) baryon number current has this property. Finally, 
this mixing is also absent if the familon has only flavor off-diagonal couplings. In all of these 
cases, the contribution of Eq. ( [4.19| ) is absent or suppressed, and the leading contributions to 



^Note that the radiatively-induced mixing operator can be written in the manifestly gauge in- 
variant form i{H^D^H — D ^H'^ H)d^^ f . This operator may be present at tree-level if the familon 
couples to "Higgs number," but we will not consider this case. 



17 



Z-f mixing come from non-logarithmically-enhanced threshold corrections which may be of 
order Such corrections are sensitive to physics at the flavor symmetry breaking scale, 
and are therefore model-dependent. 

With these caveats in mind, we now assume that the leading contribution given in 
Eq. ( [4 .191 ) is present, and determine bounds on F for various flavor symmetries. The 
logarithmically-enhanced mixing induces an effective coupling 

^cff = I A \2n^ 9zIxi9L - 9r)^ In -^dj tfj-f^ (i^Pl - sin^ Ow) ip 
2N 

= V/i \2n^ 9zIxi9L - ^r)^ In -^Ii,m^f^-f5^ , (4.20) 
[Att)^}' rriz 

where in the last step we have integrated by parts and substituted the equations of motion 
for ip. For example, for a familon coupled to tR, Nc = 3, = 1/2, and Ql — gR = —I- For 
a familon coupled to Ql, the contributions from both ti and hi must be summed. 



The effective coupling of Eq. ( 4.20 ) is constrained from various sources. For the case 
= e, a stringent constraint is provided by red giants. Familon-electron couplings lead 
to additional sources of red giant cooling, which, if too large, would destroy the agreement 
between the observed population of red giants in globular clusters and stellar evolution 
theory. Such constraints have been studied extensively in the literature [^,0 . The current 



best upper limit on the coupling is 



^ < 2.5 X 10"^^ (4.21) 

for £efj = —igfe'-)^e^ The strongest bound on the family symmetry breaking scale is for 
familon couplings that are dominantly proportional to m^, as, for example, when a familon 
is coupled only to Ir. Such results in the bound 

> 1.2 X 10^ GeV , (4.22) 

where we have taken A = F. Weaker, but still significant constraints are obtained if the 
familon coupling is dominated by rub, as when the familon couples only to Br. The bound 
in this case is 

Fi > 6.1 X 10^ GeV . (4.23) 

Notice that, in the case where the familon couples to Ir and Ql with the same charge, the 
bound of Eq. ( |4.23|) also holds for the corresponding fiavor scales. However, possibly stronger 



^For a larger coupling, the familon may be trapped in red giants and not contribute to their 
cooling. Still, they can be emitted from the Sun and change its dynamics significantly. For yet 
larger couplings, familons may be trapped in the Sun as well, but then they contribute to the 
thermal transport. Combination of these constraints exclude all couplings to electrons larger than 
this one. See Ref. [35| for further details. 
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bounds may also be possible if model-dependent non-logarithmically enhanced terms pro- 
prtional to mf are present. If the familon contribution to the induced coupling is dominantly 
proportional to m^, we find the constraint 

F^^ > 2.5 X 10^ GeV . (4.24) 

Similar bounds may be obtained from induced familon couplings to nucleons using con- 
straints from supernova SN 1987A by rescaling the bounds on axion couplings. These con- 
straints are somewhat more ambiguous because of the loss of coherence in axion emission 
due to nucleon spin fluctuations caused by scattering effects in the supernova core |^ . More 
realistic estimates were addressed in Refs. [^2i^ . The constraints yield results comparable 



to the red giant bound on electron couplings, but with larger error bars. 

Finally, we stress again that these bounds are for specific flavor symmetries. For cer- 
tain examples mentioned above, the logarithmically-enhanced contribution to the induced 
coupling is absent. For such flavor symmetries, the threshold corrections at the flavor scale 
must be studied separately for each model, and the flavor scale can be constrained only after 
the model-dependent coefficients are known. However, from the numerical estimates above, 
it is clear that the induced loop-level bounds can provide interesting constraints on flavor- 
diagonal familon couplings. Such bounds are particularly interesting for couplings to the 
third generation, as they are enhanced for large fermion masses. Note also that couplings to 
the top quark are stringently bounded and are extremely difficult to bound by other means. 



C. Bounds from effects induced by flavor-mixing 

In this section, we have so far parametrized and constrained possible familon couplings 
individually by introducing effective flavor scales, neglecting possible mismatches between 
flavor and mass eigenstates. However, as noted in Sec. ||, when bounds on a particular 
familon coupling are very stringent, such as in the case of bounds on familon couplings 
to the first generation from supernovae, one can also obtain interesting bounds on other 
familon couplings from flavor-mixing effects. In this subsection, we will consider such bounds 
in the quark sector. In extensions of the standard model with massive neutrinos, similar 
arguments hold in the leptonic sector. Additional constraints may also be obtained if the 
gauge symmetry is enlarged. 

Let us assume that the flavor and mass eigenstates coincide in the down sector. A generic 
familon coupling term for left-handed down-type quarks is then 

Cf = -^d,fdn''PLdj, (4.25) 

where / and J are generational indices. Notice that, in this section, the familon / is a real 
scalar for the diagonal couplings (/ = J), and a complex scalar for the off-diagonal ones 
(J 7^ J). (Thus, in the off-diagonal case, there is also a Hermitian conjugate term in the 
Lagrangian.) 

The up-type coupling required by gauge invariance in terms of down quark mass eigen- 
states is 
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= ^^{xud^fuj^-f5U + ■■■) , (4.26) 

where x„ = —V*jVuj- Therefore, a constraint on obtained from supernova cool- 
ing through familon-nucleon couphng imphes similar constraints on the expressions 
F^s/ i^ub^us), F,^d/ i^ub^ud) , and so on. Of course, different contributions to the same ef- 
fective coupling ^ may also have opposite signs, which must be checked in the specific 
model under investigation. 

To derive constraints on the flavor symmetry breaking scales, we must convert the quark 
level couplings to the effective nucleon-familon couplings of the form Ci^t ~ igjNN fN'-^^N . 
This can be done through a generalized Goldberger-Treiman relation. With the interaction 
given in Eq. (|4.26|) , we obtain 

9fNN = ^x.Af ) , (4.27) 
where rriN ~ 0.94 GeV is the nucleon mass, and the coefficients A^^^ are given by |43 



A(f) ~ Ai"^ ~ 0.80 , 
AW ~ aJ^ ~ -0.46 , 
A(p) ~ Ai") ~ -0.12 . 

Here, we have assumed that the flavor symmetry is anomaly-free for SU(3)c- If the fla- 
vor symmetry is anomalous under SU(3)c, there are anomaly- induced contributions to 
Eq. ( [4. 271 ); see Refs. [^,^,0] for discussions of constraints on axions. 



The effective couplings of Eq. ( [4.27|) are constrained by supernova SN 1987A. In Ref. 



the upper bound on gfpp is given as a function of gfnni Qfpp- For simplicity, we adopt the 
most conservative constraint on gjpp., 

gfpp < 3 X 10-1° , (4.28) 

and use this to estimate bounds on the flavor symmetry breaking scales. 

Under the assumption that only one down-type familon coupling exists at a time, that 

is, ignoring possible cancellations between two different contributions, we find bounds on 
third-generation couplings 

K^. > 3 X 10^ GeV ( '^""^ ( ^] , (4.29) 

'''' \3.5 X lO'V \0.80j ' ^ ^ 

F,^ > 3 X 10*^ GeV ( , ^ ( ( ^] , (4.30) 

\3.5 X 10-V VO-22 / V0-80y ' ^ ^ 

F,^, > 1 X 10^ GeV ( l^""' \ ( ( ^] . (4.31) 

'"^ \3.5 X 10-3^ \^o.98 ; vo.8oy ^ ^ 

An even stronger constraint is obtained for F^^: 
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F,^, > 8 X 10« GeV ^ ^ . (4.32) 



This constraint is, however, weaker than the laboratory bound. 

Note that the above bounds are obtained under the assumption that the mass and flavor 
eigenstates are identical in the down sector. If, on the other hand, these eigenstates were 
assumed to be identical in the up sector, familon coupling to the s- and (i-quark arises due 
to the mixing effect, and the supernova constraints on and F^^ could be used instead. 
In particular, an interesting bound is derived for F^: 



Fh>l X 10^ GeV 



(4.1 X 10-2)2 X 0.12 + (5.7 X 10-3)2 x 0.46 



(4.33) 



This bound is about one order of magnitude stronger than the bound on F^^ since, in this 
case, the effective familon coupling to the nucleon is dominated by the s-quark contribution, 
and the interaction is therefore not as highly Cabbibo-suppressed as in the case. 

Of course, there is no reason why the familon coupling is diagonalized in one sector. 
However, as the up and down sectors cannot be diagonalized in the same basis, one generally 
expects similar mixing- induced constraints in all cases. 



V. PROSPECTS FOR FUTURE PROBES: B FACTORIES 

In this section, we estimate what constraints on familon couplings may be obtained from 
the current CLEO data set and the upcoming CLEO III, BABAR and BELLE experiments. 
We make no attempt to conduct detailed experimental studies appropriate to each of these 
experimental settings. Rather, our intent here is to describe a number of analyses that are 
likely to significantly improve the present limits on familon energy scales, and, we believe, 
merit further study. We will begin with investigations of hadronic couplings of h quarks to 
familons, and then consider decays of the r lepton to familons. 



A. Bounds from B decays 

In all of the experiments mentioned above, one may search for the exclusive decays 
—>■ {Ti^,K^)f and K^f. These exclusive modes have smaller branching fractions 

than the inclusive modes b — > {d,s)f, but have clear experimental signatures due to their 

simple two-body kinematics. The form factor of 

{K-{p')\sYb{q' = 0)\B~{p)) = F,{0){p + pr , (5.1) 

which is necessary to calculate branching fractions, has been estimated by Colangelo et 
al. [^] to be Fi(0) = 0.25 ± 0.03, based on sum rules. Estimates of Fi(0) based on the 
quark model are 0.34, 0.36, 0.30, or 0.35, depending on which quark model parameters 
are assumed EH]. We could not find other estimates of this particular form factor in the 



literature, but various estimates for i? ^ tt transitions give comparable but slightly larger 
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values. This is reassuring, since they must agree in the flavor symmetric hmit. The decay 
rate B Kf is given by 



T{B Kf) 



1 m 



B „2 



(5.2) 



where /? 



m 



K 



jm 



(If the couphng is purely axial, there is, of course, no contribution 



to B ^ Kf] searches for decays to K* f are required to bound such couplings.) Neglecting 
the mass difference between the h quark and B meson, and using the naive spectator model 
for the B meson decay, one flnds 



T{B ^ Kf) 
T{h ^ sf) 



9l 



9v + 9a 



(5.3) 



The concept of the search for such exclusive decay modes is relatively simple. After 
applying the standard cuts to suppress continuum qq and lepton pair events, one looks for 
events at the T(4S') resonance that have either an isolated Kg., or an isolated charged meson 
TT^, K"^ together with large missing energy. In the center-of- momentum frame, the energy 
of the meson P = tt^, i^^, or Kg must be in the narrow range 



4 




< Ep < — 



4 



m 



mi 



m 



B , 



m 



B , 



(5.4) 



where (3 = \Jl — Am^/s = 0.0645, and mp is the mass of the meson. One can also require 
that, after excluding the isolated energetic meson whose energy is in the above range, all the 
tracks and energy deposits in the calorimeters reconstruct tub and have total energy ^/s/2 
in the center-of-momentum frame. 

Of the existing analyses, the one most similar to that described above is a search for 
5± l^vi by the CLEO Collaboration The reported upper bounds on the branching 
fractions are 1.5 x 10~^ (e), 2.1 x 10~^ (/i), and 2.2 x 10^^ (r). The reach for the Pf mode 
is expected to be worse than for ev or /iz/, as the continuum backgrounds are larger and 
the detection efficiencies are worse for mesons. However, we expect the sensitivity to such 
meson decays to be greater than to the tv mode, because the mesons have more-or-less fixed 
energy, unlike in the rv case. tt/K separation is probably difficult with the current CLEO 
data set [Q, but this analysis may still give us an upper bound on B{B iTf)+B{B — > Kf) 
somewhere at the 10~^ to 10~^ level P9| , p0| . Such a constraint would be competitive with 
the upper bound inferred from the ALEPH b svv study discussed in Sec. |III Q . 



Particle identification will be much better at CLEO III, which will allow / K separation, 
and will be even better at BABAR and BELLE. The higher luminosity at these machines 
will also help, and an upper bound of 10~^ may be possible ^9|. Such a bound would imply 
a bound on the fiavor symmetry breaking scale of F^^, F^^ ^ 2 x 10^ GeV! 



B. Bounds from r decays 

We now turn our attention to the lepton sector. The decay rate for t If is given in 
Eq. ( p.2|) . A search for t ^ If suffers from the standard model background r — Ivv. A 
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FIG. 3. (a) The e momentum spectrum from r decays. The points with error bars are our 
Monte Carlo simulation for the standard model only, normalized to the size of the current CLEO 
event sample. The solid histogram is the standard model prediction for r evu. The dashed 
histogram is the predicted spectrum with B{t ^ ef) = 3 x 10~^, again normalized to the current 
event sample, (b) The same as (a) but plotted as a ratio to the standard model prediction. 



conventional method for bounding the branching fraction to familons is to fit the momentum 
spectrum of the electron (muon) from tau decay to a linear combination of the standard 
model spectrum, which drops approximately linearly for large momenta, and a possible 
contribution from the familon mode, which is flat. The ARGUS bound quoted in Sec. [Ill A| 
was obtained by this method. CLEO has not reported a similar analysis. However, in a 
recent CLEO analysis of the Michel parameter in r decays 



511, the electron momentum 



distribution of 33531 



'z/z/)(7r 



events was presented in uniform 0.25 GeV 



bins. This distribution may be fit beautifully by the standard model alone, and contains 
about 90 events in the highest momentum bin. For reference, the contribution of a familon 
decay mode with B{t ef) = 3 x 10~^, a branching fraction near the current ARGUS limit, 
would contribute 28 events in each bin, leading to a significant excess at high momentum. 

To see the possible sensitivity given the current CLEO data set, we generated 33531 
standard model t e + missing events. The momentum spectrum from our Monte Carlo 
simulation (points with error bars), along with the predicted standard model spectrum 
(solid histogram), is shown in Fig. For comparison, we also plot the spectrum given 
a hypothetical familon branching fraction of B{t — > ef) = 3 x 10~^ normalized to the 
same number of events. In Fig. ^o, we plot the ratio to the standard model prediction 
to make the familon contribution more visible. Note that the spectrum with the familon 
contribution differs considerably in the high momentum bins. By fitting the Monte Carlo 
data to the linear combination of the standard model and familon modes, we find that 
CLEO can obtain an upper bound on the familon branching fraction of 1.6 x 10~^ (95% 
CL). This would already improve the ARGUS bound [Q on r — e/ slightly, leading to a 
lower bound of F^-e > 7 x 10® GeV. Note that in this analysis, only events with the vr^vr^z/ 
decay of the other r were used. This requirement was motivated in the original study by 
the desire to study spin correlations between decaying r pairs, but is not necessary for our 
purpose. The statistical power of our analysis may therefore be boosted by including events 
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with additional decay modes of the other r.0 

In this analysis, the systematic effects appear to be under control. The momentum 
dependence of the electron identification efficiency can be calibrated by using the actual 
data, for instance, by using radiative Bhabha events, and this calibration improves with 
statistics. In addition, the background is small in the above CLEO data sample. Indeed, 
all measurements of r decay parameters are statistically limited and we expect a similar 
situation for the familon analysis. Note also that in the above analysis we simply fit to 
the standard model contribution, allowing its normalization to vary. In principle, one can 
determine this normalization by measuring the efficiencies of r identification in each decay 
mode through methods analogous to the multi-tag methods employed in the measurement 
of Rb in Z decays We therefore conclude that a dedicated analysis could well lead to 
an upper bound on B{t ef) below the 10~^ level. 

The T ^ mode is more difficult because the muon identification efficiency is less 
well-calibrated and the statistics is slightly poorer, with 21680 events in the CLEO analysis. 
Again, however, the uncertainties are dominated by statistical errors.^ We therefore expect 
an upper bound on B{r —>■ fif) only slightly worse than that on the ef mode. 

Although we are concerned primarily with B factories in this section, we should note 
that the above analysis may also be applied at LEP. For example, in a recent OPAL analysis 
of r polarization |Q, a large sample of 25000 r ^ e events was studied. By fitting the e 
momentum distribution as described above, an upper bound on B{r ^ If) at the 2 x 10^^ 
level could be derived. (Here, we have simply scaled the CLEO results given above to the 
OPAL statistics; we have checked that the momentum distributions are sufficiently similar 
that such an approximation is valid.) Combining the four LEP experiments, we expect an 
upper bound of ~ 10~^. 

At the asymmetric B factories, the boosted center-of-momentum system would somewhat 
complicate the analysis of the electron (muon) momentum spectrum. We are not aware of 
any studies at these colliders. However, given that the event samples available at these 
machines will be much larger than the current CLEO data set, we expect BABAR and 
BELLE to place constraints significantly better than the 10~^ level, and possibly at the 
10~^ level. The study of familon decays at these machines is extremely promising, and 
worthy of further study. It must be mentioned, furthermore, that the future B factories will 
be able to improve the upper bound on the mass of Ur significantly to the level of 3 MeV, 
or possibly even 1 MeV H. As we will see in the next section, the interplay between 
bounds on the z/^ mass and bounds on branching ratios to familons is very interesting from 



^Another interesting possibility is to exploit the spin correlations between decaying r pairs by 
selecting events with u) > (see Ref. for the definition of uj). Such a selection enhances 
the right-handed t~ (left-handed t+) decaying to leptons, and thereby suppresses the electron 
momentum spectrum at the end point. The sensitivity to familon contributions at this endpoint 
is then improved. 

^In the Michel parameter analysis, a large systematic uncertainty arises because the standard 
model muon decay parameters are not assumed. For our purposes, however, we may assume the 
standard model predictions and eliminate these uncertainties. 
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a cosmological point of view. 

Finally, we note that the bounds on r — > // branching ratios are expected to be even 
better at a tau-charm factory. Ref. has shown that one can reach the level of -B(r — >■ 



ef) < 10~^ using the standard optics or even 10~^ using a monochromator. This would raise 
the lower bound on the flavor scale to Fre ^ 10^ GeV. The fif mode is more difficult, and is 
limited by the /i/vr separation capability |^5|. However, a bound better than the 10~^ level 
using a RICH detector for particle identification is expected . 



VI. IMPLICATIONS FOR NEUTRINO COSMOLOGY 

Non-standard properties of neutrinos are always interesting in cosmology, and in fact, 
heavy unstable neutrinos are advocated in certain scenarios to obtain reasonable agreement 
between theory and observation. The heavy neutrino is typically taken to be the tau neu- 
trino, and we will assume this to be the case in this section. Once a decaying neutrino 
is required, its decay into a lighter neutrino and a massless boson is the most harmless. 
Visible neutrino decays are usually severely constrained from SN 1987A. As mentioned in 
Sec. the energy released from SN 1987A was mostly carried away in neutrinos, and the 
visible luminosity of SN 1987A was much smaller. However, if neutrinos decay into visible 
particles, such as photons or electrons, neutrinos emitted from SN 1987A that decay be- 
fore reaching the earth may increase the apparent visible luminosity of SN 1987A to levels 
much larger than observed [Q. In addition, scenarios with r neutrinos decaying into 
three neutrinos are also dangerous, since, in the absence of fine-tuning, such models also 
predict large fiavor violating r decays (like r 3e) by SU(2)i gauge symmetry |T^. In 
particular, in the cosmological models to be described below, the resulting fiavor-violating 
r decay rates are already excluded by current bounds. This is because, for these three-body 
decays, the fiavor-violating r branching fraction is of order {fnr/rnv^YT~^ (to be compared 
with {fyir/my^YT'^ for the two-body familon decays), and hence is extremely enhanced for 
scenarios with neutrino masses in the currently allowed range. 

Therefore, if we adopt a massive unstable neutrino as a solution to cosmological problems, 
massless bosons are good candidates for its decay products. The familon is an example of 
such a massless boson. (In the literature, a particular example of a familon, the Majoron, 
is often considered.) If neutrinos decay to familons, cosmological scenarios, which require 
specific ranges for neutrino masses and lifetimes, then predict rates for familon signals in 
future experiments, or may even be excluded from current familon bounds, assuming an 
absence of fine-tuning. In this section, we first review some of the potentially interesting 
cosmological scenarios that require massive neutrinos. Then, assuming that the massive tau 
neutrino decays through z/,- — >• i/;/, where / = e,/x, we will discuss how well these scenarios 
may be constrained by current and future collider searches for familons. 

Among the several cosmological motivations for massive neutrinos is the "crisis" in the 
standard BBN scenario. The standard BBN scenario contains only one free parameter, 
the baryon to photon ratio r/; the abundances of the light elements are predicted once we 
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Such constraints may be evaded in scenarios with sufficiently long-lived neutrinos 



25 



fix rj. Until a few years ago, the theoretical prediction with 77 ~ 3 x 10~^° was in good 
agreement with observations. Recently, however, it has been claimed that the predictions of 
the standard BBN are disfavored by observations of the light element abundances |]58| , |59| : 
normalizing 77 with the D and '^He abundances, the observed ^He abundance is claimed to 
be smaller than the standard BBN prediction. There are several arguments against this 
viewpoint on the observational side. For example, the apparent discrepancy vanishes if one 
adopts a larger systematic error in the observed ^He abundance or if the recently 

measured D abundance in high red-shift QSO absorber systems is regarded as a primordial 
one 0.0 

On the other hand, if we regard this "crisis" as a genuine problem with the standard BBN 
theory, it can be taken as an indication of new physics beyond the standard model. There 



are several attempts to solve this crisis by a modification of the standard scenario p2| , |63| , p!6 
Here, we concentrate on a solution that uses massive unstable neutrinos to reduce the 
predicted "^He abundance. Since the *^He abundance decreases as the energy density at the 
neutron freeze out time decreases, the ^He abundance becomes smaller if N^, the "effective 
number of neutrino species" at the neutron freeze out time, is reduced. In the standard BBN, 
is 3, but it can be smaller if heavy neutrinos decay and effectively convert their energy 
density into lighter particles. For example, in the presence of the decay mode —>■ uif, if all 
the tau neutrinos are converted into thermal familons and light neutrinos, N^, ^ 2.6. BBN 



scenarios with massive neutrinos decaying to familons have been discussed in Refs. ||I3H16 
The most recent calculation shows that a massive neutrino with m,y^ ~ 10 — 20 MeV and 
T,y^ ~ 10~^ — 1 sec can resolve the conflict between theory and observation [|16|. 

Decaying massive neutrinos are also interesting for large scale structure formation. The 
standard cold dark matter (CDM) scenario, which assumes a flat universe, a scale-invariant 
initial spectrum, and that the universe is mostly fllled with slowly moving ("cold") parti- 
cles [0, is very attractive in explaining the origin of large scale structure. However, if the 



normalization of the power spectrum is fixed by the anisotropy in the temperature of the 
cosmic background radiation observed by COBE, the standard CDM scenario predicts too 
large density fiuctuations at small scales (A < 100 Mpc). Attempts to explain the scale 
dependence of the density perturbations include proposals of CDM with a small component 
of hot dark matter or with a cosmological constant pSl, or scenarios with a tilted initial 



density fiuct nation [p6 



As pointed out in Refs. |ll8|-p0|, CDM with late decaying neutrinos also provides a solu- 
tion to this problem. If the neutrino lifetime is long enough, neutrinos dominate the energy 
density of the universe at the temperature T ~ m^^. After this stage, the mass density of 
the neutrino, p^, scales as T^. Neutrinos then decay at time t ~ (py^/Afpi)^^ ~ r^^. Once 
they decay, the energy density of the neutrinos is converted to radiation energy density, 
resulting in an increase of the radiation energy density without affecting the background 
photons. This then delays the time of matter-radiation equality, the matter-dominated era 
starts later, and, with the COBE normalization, the density perturbations at small scales 
are reduced. Due to the neutrino decay, the energy density of the radiation is increased by 



^^One should note that there is another measurement of D abundances that conflicts with the one 



preferred by the standard BBN scenario [61|. Thus, this issue is still an open question. 
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the factor ~ p^{Td)/T^ ~ (m^^r;^^)^/^/Mp/^, where is the temperature just before the 
neutrino decay. As we can see, physics (approximately) depends on the combination ml^r^^. 
To obtain the correct density fluctuation at small scales, this combination must lie in the 
range {m^JkeVfiT^Jyi) ~ 50 - 150 |0 . 



In these scenarios, the neutrino mass must also lie in a specific interval. The neutrino 
must be heavier than about 50 eV; otherwise, its mass density is always smaller than (or 
at most comparable to) the mass density of the CDM, and the scenario does not work well. 
On the other hand, if the neutrino mass is above ~ 1 MeV, the neutrinos decouple from the 
thermal bath after becoming non-relativistic, and their number density is reduced. In this 
case, the constraint given above is not applicable. Furthermore, if the neutrino mass is in 
the range ~ 1 - 10 MeV while the lifetime is longer than ~ 1 sec, the neutrino mass density 
may be so large at the neutron freeze out time that ^He can be overproduced. Such lifetimes 
and large neutrino masses are therefore also disfavored from BBN considerations [^. In 
this section, we consider the mass range 50 eV ^ m,^^ < 10 MeV, with the above caveats in 
mind. 

To summarize, we consider the following two cosmological scenarios: 

• BBN : m^^ ~ 10 - 20 MeV, and r^, ~ 10"^ - 1 sec. 

• Structure formation : (mi,^/keV)^(r,^^/yr) ~ 50 — 150, with 50 eV ^ m^^ < 10 MeV. 

These scenarios require decays to familons ^ ^if- As we discussed previously, this 
process is related to the decay modes r ^ If and b ^ sf though SU(2)j, and GUT gauge 
symmetries, respectively. Thus, searches for these decay modes are interesting tests of these 
scenarios. 

Let us start with the r familon decay mode. If the decay mode Ur — ^ ^if exists, by SU(2) l 
symmetry, the charged r lepton must also have flavor-changing couplings to familons: 

Cf = jd,f {gr^yrYPm + gl'fYPLl) + h.c. (6.1) 

If there is no fine-tuning, gjj ^ g'l'^'- Notice that, if the right-handed leptons also transform 
under the flavor group, they also couple to familons, and such interactions may increase 
the rare r decay rate. The following argument is therefore conservative. From the above 
Lagrangian, we obtain the decay rate 

Tfl 

T(t If) = , (6.2) 

and, using Eq. ( [4.14| ), we find 

B{t ^ //) = 

^8.1xl0-x (^-^1 (^M^l . (6.3) 




In Fig. ^ contours of constant B{t If) are shown in the {m^^,T^^) plane, assuming 
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FIG. 4. Contours of constant B(t — > //) (10^'^, 10~^, and 10^'', from below), assuming 
^utui — natural 811(2)/, gauge relation in the absence of fine-tuning. The regions pre- 

ferred by the BBN and structure formation scenarios discussed in the text are also shown. Note 
that the parameter region for structure formation with rriy^ > 1 MeV may be unreliable because 
of limitations in the approximations used to derive the preferred region. The lightly shaded region 
is disfavored by BBN. 



From Eq. ( |6.3D , we see that, in the absence of fine-tuning, the BBN scenario with decaying 
neutrinos predicts B{t If) ~ 10"'^ to 10"''. The current bounds B(t~ /^~/) < 
4.6 X 10-3 (95% CL) and B{t- f) < 2.6 x lO^^ (95% CL) § therefore do not constrain 
this scenario. However, if the sensitivity of future experiments is improved by one order of 
magnitude or more, the predictions of this scenario may be tested, and, if the scenario is 
correct, exotic r decays may be seen. (See Fig. ^.) 

The scenarios motivated by structure formation are also interesting. In this case, 
Eq. ( |6.3[ ) implies B{t —>■ If) ~ 10"^ to 10~^. Comparing this result with the current 
bound, part of the parameter region of this scenario is already excluded. As discussed in 
Sec. |V B| , future experiments may reach a sensitivity for B{t —>■ If) of 10"^ or possibly 10~^. 
Thus, if the CDM scenario were realized, the familon decay mode is likely to be found if the 
neutrino is lighter than ~ 1 — 10 keV. The region disfavored by BBN jl^ is also shown by 
a light shading. 

Up to now, we have only used SU(2)i gauge symmetry to relate the neutrino-familon 
interaction to existing and future constraints on r decays. However, if we assume that the 
same familon also couples to down-type quarks, the above cosmological scenarios may also 
be probed by rare b decays. Such is the case in SU(5) GUTs, where the lepton doublet and 
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FIG. 5. As in Fig. but with contours of constant B{b ^ sf) (10~^, 10~^, and 10^^, from 
below), assuming F^^^^ ~ F^, as would be the case in GUTs without fine-tuning. 



right-handed down-type quarks are in the same multiplet, and so we also have a couphng of 
the form 

C D ^d.fg'^h'PRS + h.c. , (6.4) 
and similarly for d. With this Lagrangian, we obtain 



B{h sf) 




2 . . / 1 MeV\^ (\ sec 




(6.5) 



where we have used = 4.5 GeV, and = = 1.6 x 10^^^ sec. Comparing Eq. 
with Eq. ( |6.3| ), we can see that the branching ratio B{b ^ sf) is enhanced by about two 
orders of magnitude relative to B{t If). This results from an enhancement by a factor 
(nib/mr)^, and also the fact that the total decay rate of the b quark is Vcb suppressed, and 
so is even smaller than that of the r lepton, despite its larger mass. Contours of constant 
B{b —>■ sf) are shown in Fig. |^. For the cosmologically-motivated scenarios, the ranges of 
B{b ^ sf) are 10^^ to 10^^ (BBN), and 1 to 10^^ (structure formation). Because the b 
decay rates are so enhanced, in each case, the high value in the predicted range is above 
bounds which can be expected from the current data. Future experiments may reach a 
sensitivity of -8(6 — > sf) ~ 10~^ (corresponding to B{B ^ Kf) ~ 10~^ in Sec. |V A| ) and 
will thus cover most of the preferred parameter regions. 
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Finally we comment on m^^ measurements at future B factories. As we mentioned earlier, 
the upper bound on m^^ will be significantly improved at future B factories, to the level 
of 1 - 3 MeV. The parameter region that will be covered has a significant overlap with the 
neutrino mass required in the above mentioned cosmological scenarios, and hence such mass 
measurements provide another probe of these scenarios. The BBN scenario with massive 
unstable neutrinos will be fully tested by the tau neutrino mass measurement in future B 
factories. On the other hand, for the structure formation scenario, most of the interesting 
parameter region [m,^^ < 1 MeV) may be covered by the search for b —>■ sf, and even if the 
neutrino mass is above ~ 1 MeV, this scenario can be checked by the direct measurement 
of the tau neutrino mass (though this region is disfavored by BBN). Therefore, in this case, 
measurements of the mass and the branching ratios will have complementary roles. 

VII. CONCLUSIONS 

If global family symmetries play a role in determining the patterns of masses and mixings 
of the quarks and leptons, they must be spontaneously broken. Familons (or Majorons), the 
massless Goldstone bosons associated with these broken symmetries, allow rare opportunities 
to probe the physics at very high mass scales in a multitude of low energy settings, and their 
discovery will signal a breakthrough in attempts to understand the flavor structure of the 
standard model. 

The experimental investigation of familons has in the past focussed on familons coupled 
to the first two generations. As we reviewed, such investigations have led to stringent lower 
bounds on the flavor breaking scale F of ~ 10^ GeV and ~ 10^^ GeV in the leptonic 
and hadronic sectors, respectively. In contrast, bounds for familons coupling to the third 
generation are much less thoroughly studied. In the lepton sector, constraints from rare 
r decays lead to constraints F ^ 10^ GeV; in the hadronic sector, no bounds have been 
previously reported. The lack of study of third generation familons is conspicuous, especially 
in light of their cosmological relevance and the upcoming B factory experiments, which hold 
promise for studying b and r decays with great precision. 

Motivated by these considerations, we have presented a large and eclectic group of bounds 
which we believe place the most stringent constraints on third generation familon couplings. 
As emphasized in Sec. |T1|, the experimental and astrophysical implications of familons vary 
strongly with the underlying flavor symmetry and depend on whether the flavor symmetry 
is real or complex, the familon couplings are axial or vector-like, and whether they are 
flavor- diagonal or non-diagonal. For instance, bounds on K decay in the presence of mixing 
effects may in general impose bounds of order 10^ GeV on third-generation flavor breaking 
scales, but there are classes of models in which such bounds do not apply. It is therefore 
important to consider a wide variety of experimental signatures. Different signatures are 
also related through some well-motivated theoretical considerations: in the absence of fine- 
tuning, the familon couplings of particles related by gauge symmetry are expected to be 
similar in strength. Probes of r decays to familons are therefore indirect probes of i^r 
familon couplings, and in SU(5) grand unified theories, these are both related to b decays. 

We began by considering bounds from currently available accelerator data. Present values 
for neutral meson mass splittings imply bounds on the flavor scale of ~ 10^ to 10^ GeV for 
real familons. Bounds from neutral meson decays to leptons are significantly weaker, at the 
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level of 10^ GeV, and require both hadronic and leptonic couplings. More promising are 
bounds from exotic b decays at LEP. By extrapolating from current bounds on 6 — >■ sh'P, 
we estimate that an analysis of currently available LEP data may provide a sensitivity to 
B{b sf) at the level of 1.8 x 10"'^, leading to probes of flavor scales of the order of 10^ 
GeV. 

Familons also have astrophysical implications, as they may lead to anomalously fast 
cooling of supernovae, red giants, and white dwarfs. Bounds from direct couplings to Ur 
are generally weak. However, couplings of familons to particles of the third generation may 
also induce couplings to electrons and nucleons radiatively or through flavor mixing effects. 
Bounds on such couplings are model-dependent, but may be stringent; in the simple case 
where a familon couples diagonally to t quarks, a bound of F > 10^ GeV from radiatively 
induced couplings may be set. 

Finally, having evaluated a host of new bounds, we considered the prospects for analyses 
at future B factories. Such colliders are ideal experimental environments for searches for 
rare r and b decay modes and are expected to have greatly improved statistics. We find that 
probes of branching fractions of 10"'^ (10^^) for r (6) decays may be possible. As discussed 



in Sec. under the assumption that the flavor scales for these decays are naturally related 
to the scales for z/,- couplings, such precise probes are sensitive to parameter regions favored 
by various BBN and structure formation scenarios, where a massive unstable neutrino is 
motivated by possible discrepancies in cosmo logical data. In fact, parts of the parameter 
regions in such scenarios are already excluded, and future searches will be able to explore 
large portions of the cosmologically-favored parameter space. Given the present lack of 
analyses studying third generation familon couplings, studies at all of these experiments, 
and particularly the B factories, are strongly encouraged. 

Note added. After the completion of this work, Ref. appeared, in which the impact of 
decays Ur —>■ Vef on BBN was examined. Tau neutrino masses of 0.1 to 1 MeV and lifetimes 
of 50 to 2 X 10'^ s were found to be allowed, implying rates for r — > e/ and b ^ df within 
reach of future experiments. (See Figs. H and ||.) 
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